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ABSTRACT 

The Galaxy Evolution Explorer (GALEX) satellite has obtained high time resolution ultraviolet 
photometry during a large flare on the M4 dwarf star GJ 3685A. Simultaneous NUV (1750 - 2800 A) 
and FUV (1350 - 1750 A) time-tagged photometry with time resolution better than 0.1 s shows that 
the overall brightness in the FUV band increased by a factor of 1000 in 200 s. Under the assumption 
that the NUV emission is mostly due to a stellar continuum, and that the FUV flux is shared equally 
between emission lines and continuum, then there is evidence for two distinct flare components for 
this event. The first flare type is characterized by an exponential increase in flux with little or no 
increase in temperature. The other involves rapid increases in both temperature and flux. While the 
decay time for the first flare component may be several hours, the second flare event decayed over less 
than 1 minute, suggesting that there was little or no confinement of the heated plasma. 
Subject headings: stars: flare — stars: late-type — stars: individual (GJ 3685A) — stars: variables: 
other (UV Ceti) — ultraviolet: stars 



1. INTRODUCTION 

UV Ceti stars exhibit many of the activity phenom- 
ena observed on the Sun, such as flares, dark spots 
and variable emission fr om the chromosphere and corona 
l)Thomas fc Weissj l2004). Such surface activity is linked 
to their magnetic field strength and surface coverage, 
which can be several times larger than that found on 
the Sun. When quiescent, UV Ceti flare stars are low 
luminosity (My > 8), low temperature (T = 2500 - 
4000 K) M-type dwarfs. During outbursts, these stars 
brighten dramatically over timescales of a few sec- 
onds to hours at all wavelength s from X-ray to radio 
(Hai sch. Strong fc Rodonclll991J) . 

The vast majority of stellar flare observations involve 
photometric moni toring using the Johnson UBV filters 
l)Lacv et al.lll97tt) . Normally this involves high speed 
sampling using a single transmission filter. Occasionally, 
however, observations are obtained in multiple filters by 
sequencing using time steps of a few seconds to about 1 
minute. This technique has the advantage of obtaining 
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flare colors, allowing estimates of temperatures and sizes, 
and is reasonably valid during the long duration decay 
phase of the events. However, during the impulsive phase 
the fluxes can change on time scales of seconds, so the 
color data becomes highly uncertain. A second limitation 
is that the flare energ y distribution peaks at wa velengths 
shorter than 4000 A |van den Oord et al.lll99^ . so that 
the Johnson photometry only samples the wing of this 
distribution, which is contaminated by the stellar photo- 
sphere. 

The GALEX satellite l)Martin et all2005|) provides the 
opportunity for obtaining vastly improved photometric 
measurements of stellar flares. Using a dichroic beam 
splitter and two photon counting detectors, this telescope 
is able to simultaneously monitor a 1.25° field of view 
in the FUV (1350 - 1750 A) and NUV (1750 - 2800A) 
regions with 5 arcsec and 6.5 arcsec angular resolution 
respectively. This allows the accurate determination of 
UV fluxes and colors at a time resolution which is lim- 
ited only by the count rate and the required S/N values. 
In this paper we describe a serendipitous observation of 
a large flare seen by the GALEX telescope on the dM4e 
star GJ 3685A. Some of the physical properties of this 
event are discussed and we develop a schematic model 
which is compared with models developed from other 
flare observations. 

2. OBSERVATIONS AND DATA REDUCTION 

The flare event on GJ 3685A was detected on 2004 
April 20 during a simultaneous FUV and NUV ob- 
servation of the Medium Imaging Survey (MIS) field 
MISDR1_13062_0283. The exposure began at 19:42:06 
UT and lasted for 1244s. At present, the standard 
GALEX Data Analysis Pipeline is only designed to pro- 
duce a calibrated image of the field and thus does not 
retain timing information. To construct a photometric 
time sequence, we began with the original time-tagged 
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photon lists and integrated all photons within a 1.8 ar- 
cmin aperture centered on the source. The background 
count rate was measured within an annulus extending 
from 1.8 to 2.4 arcmin around the star, during the first 
500 s of the observation, in order to avoid contamination 
by scattered light when the flare was at its brightest. The 
FUV and NUV background rates are 3.2 cts s^ 1 and 38.5 
cts s _1 , respectively. 

The G^LiSV microchannel plate detectors exhibit a lo- 
cal non-linearity of 10% affecting point sources with in- 
put co unt rates of 90 cts s - 1 in FUV and 470 cts s _1 in 
NUV l|Morrissev et al.ll2005D . In this paper, we use Mor- 
rissey's published local non-linearity corrections, which 
have been calibrated with white dwarf standard stars up 
to 3300 cts s" 1 in FUV and 8400 cts s" 1 in NUV. The GJ 
3685A flare exceeded this calibrated range in the FUV 
for a few seconds during the brightest parts of the flare: 
we have omitted these data from the analysis because 
the uncertainty in the true count rate is unknown in this 
high count-rate regime. 

The results of the analysis are given in Figure ^ which 
shows both the flux variations and the ratio of the FUV 
to NUV fluxes. In these plots the photon count rates have 
been converted to fluxes (in erg cm -2 s _1 A -1 ) using 
the conversion factors of 1.4 x 10 ~ 15 and 2.6 x 10~ 16 fo r 
the FUV and NUV, respectively l|Morrissev et alj f2005l. 
Initially, the star showed a 'quiescent' flux of ~ 10~ 15 erg 
cm~ 2 s _1 A -1 and a flux ratio (or color) of 1. The start 
of the stellar flare is lost in the noise of the background. 
A reasonable estimate seems to put the start somewhere 
between 450 and 500 s and it is definitely underway in 
both the NUV and FUV bands by 500 s. Between 500 s 
and 630 s the flux in both bands increases exponentially, 
with an e-folding time of ^60 s. During this time the 
average color remains constant, at a value of near 0.8. 

At 650 s the rate of flux increase jumps dramatically, 
with e-folding times of 23 s in the NUV and 12 s in the 
FUV, and the FUV/NUV ratio increases from < 1 to 
more than 6. Both temperature and FUV/NUV ratio 
reach a peak at 700 s. During this peak period, we have 
omitted the data because the FUV non-linearity correc- 
tion is uncertain. 

After the peak, which only lasts for about 20 s, the 
flux and FUV/NUV ratio rapidly decrease. Note that 
the time variations of both the color and the flux for the 
100 s centered at the flare peak are symmetrical about 
the peak, i.e. the rise times and decay times are about 
the same. By 750 s the color ratio is again near 1 and 
the flux settles into a slow decline which is interrupted 
at 865 s by a second flux enhancement. This second en- 
hancement is somewhat more complicated than the first, 
but follows the same basic trend of rapid flux increase ac- 
companied by a rapid increase in FUV/NUV, this time 
to values in excess of 10. The decay phase for the second 
enhancement is slower and more complex than was the 
case for the first enhancement. By the end of the obser- 
vation the flux and color had returned to the levels seen 
prior to the second enhancement, though the flux is still 
well above the 'quiescent' level. 

3. ANALYSIS 

3.1. Pre-flare Activity 

On 2004 March 2 GALEX obtained a 110 sec 'qui- 
escent' observation of GJ 3685A which showed a flux of 



1.7 ± 0.5 x 10" 16 erg cm" 2 s" 1 A^ 1 in the FUV and 
3.0 ± 0.3 x 10~ 16 erg cm~ 2 s" 1 A" 1 in the NUV, giv- 
ing a color ratio of 0.57 ± 0.19. On 20 April, before 
the onset of the flare event the average 'quiescent' flux 
was 1.54 x 10~ 15 erg cm" 2 s _1 A -1 in the NUV and 
1.52 x 10~ 15 erg cm -2 s _1 A -1 in the FUV band, result- 
ing in a flux ratio of 1. Since the FUV/NUV flux ratio 
is directly related to the effective emission temperature 
(see section l3~2*| we see that preceding the flare the stel- 
lar emission had not only increased in flux by a factor of 
5 in the NUV, but also increased in temperature. 

One possible explanation is that the pre-flare enhance- 
ment results from an increased level of microflaring ac- 
tivity. To test this notion, we performed an analysis on 
the first 400 s of t he data set using the binnin g tech- 
nique described by lE.obinson et alJ l) 19951 11999ft . Indi- 
vidual microflares can be detected by binning the time 
sequence and searching for times when the counts per bin 
are significantly higher than those expected from random 
chance. The maximum visibility occurs when the binning 
factor is near the lifetime of the event. 

A distribution of small events which are not individ- 
ually time resolvable can be determined by comparing 
the distribution of count rates with a reference distribu- 
tion from a non- variable source, as described in Robinson 
(1999). This method reveals that when the binning fac- 
tor is much smaller than the flare lifetime, then event 
noise will dominate and both distributions will be Pois- 
son. However, as the binning factor approaches the flare 
lifetime, the microflares show up as an enhancement of 
the high count rate tail of the distribution. 

Examining the time sequences for binning factors rang- 
ing from 1 s to 20 s for the GJ 3685A flare data shows 
no strong evidence for either individual or unresolved 
flare events, while a regression analysis showed that the 
level of activity remained statistically constant prior to 
flare onset. A problem, however, is that the background 
levels are large compared to the stellar signal. For ex- 
ample, in the NUV the estimated background is 38.5 
counts s _1 , while the stellar signal was only 7.5 counts 
s _1 . Thus, noise from the background would swamp all 
but the largest microflare events. 

3.2. Empirical Modeling 

Since we have no spectroscopic information about 
this flare it is not possible to perform any type of 
detailed modeling or radiative transfer calculations on 
these UV data. However, some good qualitative and 
semi-quantitative insights can be obtained by assuming 
black-body emission. The reader should be aware that 
this assumed model is only one of several possible flare 
model scenarios. However, this simple theoretical ap- 
proach has been widely used in the interpretation of op- 
tical photometry of flare events and yielded some impor- 
tant insights into flare phenomena llKahler et all 119821 
Ivan den Oord et al.lll996t |Hawlev et al.ll2003|) . We have 
a definite advantage in these UV observations since the 
photospheric emission from the star is so small that it 
can be neglected. 

The first step in our simplified analysis is to determine 
the relation between the black body temperature and 
the measured FUV/NUV flux ratio. We have calculated 
the black body spectrum for a number of different tem- 
peratures, multiplied by the effective area curves for the 
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GALEX FUV and NUV filters and then divided by the 
total effective area for each filter. This gives the expected 
flux for a black body, which can be directly related to the 
average flux deduced from the measured count rates. The 
relation between temperature and flux ratio is shown in 
Figure El Note that the flux ratio increases dramatically 
between 4000 and 19000 K (when the peak of the black 
body distribution falls within the FUV bandpass) and 
then increases much more slowly toward high tempera- 
tures. The ratio eventually reaches a value of just over 
3 at temperatures of 50,000 K. This immediately shows 
that our black-body assumption is not totally valid, since 
we measured ratios during the flare peaks in excess of 6. 
This is not a surprise, since the FUV wavelength region 
has numerous strong emission lines (e.g. C II, C IV, Si II, 
Si IV, etc) which can dominate the continuum, especially 
during 'quiescent' periods outside of major flare events 
l)van den Oord et al]ll996|) . 

Assessing the relative contributions from continuum 
and/or emission lines at near U V wavelengths is still a n 
open problem for flare research (Ha wlev fc Fisherl ll992'). 
We note that the GALEX NUV channel is not signifi- 
cantly contaminated by the (expected) strong Mgll emis- 
sion lines at 2800A due to the low transmission of the 
GALEX NUV instrument at this wavelength. However, 
previous flares have shown an (albeit far less) increase 
in the NUV emission fro m the Fe ll line multiplet at 
2600A llHaisch et al.lll987j) . whereas l&utler et all l)1981fl 
have observed a significant rise in the near UV contin- 
uum during a flare on the star Gl 867A. For the present 
case of the flare on GJ 3685A, it is difficult to imagine 
how line emission could be the sole major contributor 
to the more than 5 UV magnitude increase in NUV flux 
observed during the flare. Thus, in determining the effec- 
tive temperature from the FUV/NUV flux ratio, we ar- 
bitrarily assume that approximately half of the FUV flux 
was from emission lines and the other half from black- 
body emission. This is consistent wit h the time-resolved 
spectroscopy of lHawlev et alJ l)2003D , where the FUV 
line-to-continuum flux ratio is 0.8, during the brightest 
flare on AD Leo report ed in that paper. Furthermore, 
lHawlev fc Fisherl l)1992|) also provide support for the no- 
tion that the GALEX NUV band is dominated by the 
star's continuum. 

A black body of given temperature has a well de- 
fined emission. Thus, having determined a temper- 
ature we can define the effective stellar surface cov- 
erage required to produce the observed NUV flux at 
the earth, assuming a distance to the star of 14.6 pc 
ijAstronomisches Rechen-InstitutJ Il998l) . The results of 
the calculations are presented in FigureOl where we have 
converted the flare surface area to an effective source ra- 
dius. The calculated temperature in the pre-flare phase 
was about 12,000 K with a deduced scale of about 3500 
km. The exact values are uncertain because this phase is 
dominated by the emission lines in the FUV and the ac- 
tual structure is most likely composed of small elements 
spread over a large area. When the flare starts the ef- 
fective temperature drops to about 10,000 K and the 
radius increases linearly from 3,500 km at 500 s to about 
25,000 km at 650 s, suggesting an initial expansion ve- 
locity of at least 140 km s _1 . The deduced temperature 
increases dramatically at 660 s, becomes undefined near 



the flare peak and then rapidly decreases back to a value 
of about 10,000 K. During this time the deduced radius 
decreases to 7,000 km and then increases back to 25,000- 
28,000 km, about the same value as that seen before the 
large temperature increase. The temperature and area 
then remain relatively constant until the start of the sec- 
ond temperature enhancement at 880 s, where the area 
again dramatically decreases and then slowly recovers to 
a value of about 32,000 km. 

Having estimated a temperature and area for the flare 
it is possible to determine the total optical and UV en- 
ergy by integrating under the appropriate black body 
curve. The result of this calculation is shown in the Fig- 
ure 0] The overall appearance of this event is two short 
duration bursts at 670-730 s and 880-970 s superimposed 
on a more slowly varying event which started at 450 s and 
which was slowly decaying at the end of the observation. 
The flare is obviously very energetic, with luminosities 
reaching more than 10 32 erg s _1 and a total integrated 
energy in excess of 10 34 erg. This is comparable to the 
largest recorded M-dwarf flare ou tbursts, such as the 
1985 flare on the star AD Leonis ijHawlev &: Pettersonl 
I1991|) . Another way to put this into perspective is to 
note that the total (optical plus UV) luminosity emit- 
ted from the photosphere of a M4 star is on the order 
of 2.0 x 10 31 erg s _1 . Thus, at the peak the flare would 
outshine the star by more than a factor of 10. 

4. DISCUSSION 

Assuming that our black body emission model dis- 
cussed in Section 3.2 provides a reasonable description 
of the flare event on GJ 3685A (and given the caveat 
that there may be other models that can provide an al- 
ternate interpretation of the data that we have chosen 
not to explore due to the lack of spectroscopic informa- 
tion on this flare), we can now postulate a schematic 
model for the event which satisfies all of the observa- 
tions and is consistent with observations of flares seen 
on other stars. The flare itself appears to be composed 
of two distinct types of events, which we will refer to 
as type A and type B. The flare begins with a type 
A event, which is characterized by a FUV/NUV flux 
ratio that is nearly the same as that seen in the pre- 
flare activity (plage) and remains constant as the flux 
increases. We conclude that the spectrum has a very 
plage-like appearance and that the source structure, and 
possibly the heating processes, are also very plage-like. 
We postulate that this phase was initiated by a mag- 
netic reconnection somewhere within a highly stressed 
active region. This initial energy release triggers recon- 
nection events i n adjacent magnetic structures, leading 
to an avalanche l)Lu fc Hamiltonll9911:lCarbonneau et alJ 
|2001[) which propagates throughout the region at a veloc- 
ity of at least 140 km s _1 , which probably represents the 
local Alfven velocity. The increase in flux is then simply 
the result of expanding area coverage. 

At about 670 s (roughly 3 minutes after the start 
of the flare) the disturbance which is responsible for 
flare A intersects a highly unstable magnetic structure 
and triggers an explosive release of energy, resulting in 
the first type B flare (Bl), which occurs between 670 
s and about 730 s. The rapid increase in temperature 
and flux is compatible with a chromosphe ric evapora- 
tion model l)Haisch. Strong. &: Rodonolll991j) . in which a 
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large amount of energy (probably in the form of energetic 
particles) is suddenly released near the top of a magnetic 
loop complex and propagates down towards the photo- 
sphere. When it strikes the denser atmospheric layers 
it impulsively heats the material, which then expands 
back into the loop. Normally, this type of event will 
have a long lasting decay as the heated material cools 
through conduction and radiation. In this case, however, 
the short duration of the decay suggests that the energy 
release wa s sufficient to cause a disruption of the confin- 
ing loops (Realc, Bocchino. &: Pereal2002j) . The source 
size for flare Bl is also significantly smaller than that 
of flare A, since the deduced flare area decreases rapidly 
when the emission from Bl dominates that from flare A. 

The structure responsible for flare Bl is probably near 
the edge of the active region associated with flare A, since 
the expansion of flare A ends during flare Bl. Appar- 
ently, all of the available structures which can support 
flare A emission have been activated by that time. It 
is possible that the disturbance responsible for flare A 
could continue into regions outside of the initial active 
region. More likely, the energy release initiating flare Bl 
also generated an Alfv enic shock, similar to t he Morton 
waves seen on the Sun ijAthav fc MortonH 961). This dis- 
turbance propagates to an unstable magnetic structure 
in a nearby active region and triggers a second B type 
flare (flare B2) which starts at around 880 s. This flare 
is both stronger and more complex than Bl, indicating 
more complex field configuration. It also shows a better 
developed thermal decay, implying that at least some of 
the hot plasma is confined in the loops. Note also that 
this is an isolated event, without an accompanying type 
A event, since there was no increase in source size prior 
to the onset of the temperature enhancement. 

By the end of the observation it appears that flare B2 
has completely faded and flare A is again visible. Un- 
fortunately, we do not have much information about the 
decay phase of flare A. However, from the fact that the 
integrated flux at the end of the observation at 1250 s is 
very similar to that at 800 s, after the decay of flare Bl, 
suggests that it may take an hour or more for flare A to 
decay to pre-flare levels. 

5. CONCLUSIONS 

We report on a large flare from the dM4e star 
GJ 3685A which was simultaneously observed in the 
NUV and FUV by the GALEX satellite. Under the as- 



Astronomisches Rechen-Institut 1998, ARI Database for Nearby 

Stars, http:/ /www. ari.uni-heidclbcrg.de/aricns 
Athay, K.G., & Morton, G.R 19&1, ViS, 935 
Butler, C.J., et al. 1981, MNRAS, 197, 815 

Carbonneau, P., Mcintosh, S.W., Liu, H-L, & Bogdan, T.J. 2001, 

Solar Physics, 203, 321 
Gliese, W., & Jahreiss, H. 1991, Preliminary Version of the Third 

Catalogue of Nearby Stars. Greenbelt, MD: Goddard Space 

Flight Center 
Haisch, B., et al. 1987, A&A, 181, 96 

Haisch, B., Strong, K., & Rodono, M. 1991, ARA&A, 29, 275 
Hawley, S.L., & Fisher, G.H. 1992, ApJS, 78, 565 
Hawlcy, S.L., & Petterson, B.R. 1991, ApJ, 378, 725 
Hawley, S.L., et al. 2003, ApJ, 597, 535 
Kahler, S., et al., ApJ, 252, 239 

Lacy, C.H., Moffett, T.J., & Evans, D.S. 1976, ApJS, 30, 85 



sumption of a blackbody emission model in which half 
of the FUV flux arises from line emission and half from 
continuum, and that the NUV flux arises solely from an 
increase in the continuum, there is strong evidence for 
two distinct classes of flares during this event. The first 
(type A) is characterized by an exponential increase in 
flux with little of no change in emission temperature, as 
measured from the FUV/NUV flux ratio. This behavior 
is compatible with an avalanche model. The fact that 
the flux ratio during this phase was very similar to that 
seen in the pre-flare plage emission also supports the idea 
tha t plages are h eated by microflares, as first proposed 
by ilPnrkerllli)*,^,. 

The second class (type B) is characterized by an impul- 
sive increase in both flux and temperature. This class is 
consistent with a classical explosive event and may be re- 
lated to the solar two ribbon flare. Two flares of this type 
were seen during the event. The first (flare Bl) probably 
originated in the same active region as the type A flare. 
The lack of a well developed decay phase suggests that 
the confining magnetic structure was disrupted. The sec- 
ond flare (B2) may have originated outside of that region, 
being triggered by a disturbance formed during the flare 
Bl energy release. 

The 2004 April 20 flare began when GJ 3685A was al- 
ready in a state of enhanced activity in which both the 
flux and temperature are substantially higher than in a 
previous GALEX observation on 2004 March 2. It is un- 
clear whether this enhancement is the result of coronal 
variations, or arises from the rotation of an intense active 
region onto the visible surface of the star, or if it comes 
from some pre-flare energy release. In order to obtain 
better models of future dMe star flare events we recom- 
mend using the low-resolution spectroscopic grism mode 
of GALEX, instead of the imaging photometric mode. 
Such future observations would provide better insights 
into the relative contributions from the UV continuum 
and/or line emission during these large releases of stellar 
energy. 
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Fig. 1. — (Top) Ratio of FUV to NUV fluxes. In this and all other figures, data are omitted during the time intervals in which the 
FUV count rate exceeds the maximum calibrated count rate of 3300 cts s _1 . (Bottom) Calibrated fluxes from the GALEX FUV (solid) 
and NUV (dashed) channels for the entire observation of 2004 April 20. Data have been binned into 25 s intervals for 0-450 s, 5 s intervals 
between 450 and 630 s and 1 s intervals for times greater than 630 s. 




Fig. 2. — Relation between black body temperature and the FUV/NUV flux ratio determined by convolving the appropriate black body 
distribution with the GALEX FUV and NUV filter effective area curves. 
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Fig. 3. — (a) Effective black body temperatures derived from the measured FUV/NUV flux ratio using the calibration presented in 
Figure|21 Preflare activity has been estimated using an integration over the first 400 s. Data between 400 s and 630 s is deduced using 5 s 
bins, while later data uses 1 s bins, (b) Effective radius of the source deduced from the observed flux, as described in the text. Data for 
the first 630 s has been binned into 5 s intervals, later data has been binned into 1 s increments, (c) Measured FUV flux, binning as in (b). 
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Fig. 4. — Estimated luminosity, assuming black body emission at the temperature and source size presented in Figure l3l 



